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Density functional studies of 26K superconducting LaFeAs(0,F) are reported. We find a low 
carrier density, high density of states, N(Ef) and modest phonon frequencies relative to T c . The 
high N{Ef) leads to proximity to itinerant magnetism, with competing ferromagnetic and antiferro- 
magnetic fluctuations and the balance between these controlled by doping level. Thus LaFeAs(0,F) 
is in a unique class of high T c superconductors: high N(Ep) ionic metals near magnetism. 

PACS numbers: 74.25. Jb,74.25.Kc,74.70.Dd,71.18.+y 



Understanding the interplay between superconductiv- 
ity and spin fluctuations especially near magnetic phases 
is an important challenge. The discovery of a new family 
of layered superconductors containing the magnetic ele- 
ments Fe and Ni and critical temperatures up to T C =26K, 
specifically LaOFeP (T C =4K, and ~7K with F doping)^ 
LaONiP (T C =3K), 3 and F doped LaOFeAs with T C =26K 
(Ref. i) raises questions about the relationship between 
magnetism and superconductivity, the origin of the re- 
markably high T c , and the chemical and structural pa- 
rameters that can be used to tune the properties. Here 
we show that LaOFeAs is in fact close to magnetism, with 
competing ferromagnetic and antiferromagnetic fluctua- 
tions, with the balance controlled by doping. Further- 
more, we identify La(0,F)FeAs as a nearly 2D, low car- 
rier density metal, with modest phonon frequencies rela- 
tive to T c , and high density of states. 

The crystal structure is layered with apparently dis- 
tinct LaO and transition metal pnictide layers, (see e.g. 
Ref. 0) Importantly, it forms with a wide range of rare 
earths and pnictogens, with the transition elements Mn, 
Fe, Co, Ni, Cu and Ru, and additionally, related fluorides 
and chalcogenides are knowni^^ It may be that there 
are other superconducting compositions that remain to 
be discovered in this family. 

Our calculations were done using the local spin den- 
sity approximation (L(S)DA) and the generalized gradi- 
ent approximation (GGA) of Perdew, Burke and Ernzer- 
hof (PBE)£ We used the linearized augmented planewave 
(LAPW) method with tested basis sets and zone sam- 
plings for the electronic structured Two independent 
LAPW codes were used - an in-house code, and the 
WIEN2K package. 9 Transport properties were calculated 
using BoltzTraP^ The phonon dispersions were ob- 
tained in linear response via the quantum espresso code 
and ultrasoft pseudopotentials with a cutoff energy of 50 
Ry and the PBE GGA. The codes were cross-checked by 
comparing band structures and zone center Raman fre- 
quencies. We used the experimental lattice parameters 
a=4.03552A, c=8.7393A, for LaFeAsO. The internal co- 
ordinates were determined by LDA energy minimization 
as ZLa=0.1418 and za s =0.6326. The corresponding Ra- 
man modes have mixed character and frequencies of 185 



cm -1 and 205 cm -1 . The corresponding ultrasoft pseu- 
dopotential values were 181 cm -1 and 203 cm -1 . The 
FeAs layers consist of a square lattice sheet of Fe coor- 
dinated by As above and below the plane to form face 
sharing FeAs4 tetrahedra. These are squeezed along c 
(the As-Fe-As angles are 120.2° and 104.4°). The Fe-As 
distance is 2.327A. The Fe-Fe distance is 2.854A, which 
is short enough for direct Fe - Fe hopping to be impor- 
tant, while the As - As distances are 3.677A, across the 
Fe layer, and a=4.036A, in plane. 

The electronic density of states (DOS) and band struc- 
ture are shown in Figs. [1] and [2J The band structure is 
similar to that reported for LaFePO^ii Near the Fermi 
energy, Ep it is well described as coming from 2D metal- 
lic sheets of Fe 2+ ions in an ionic matrix formed by the 
other atoms. There is a group of 12 bands between -5.5 
eV (relative to Ep) and -2.1 eV. These come from O p 
and As p states, with the As p contribution concentrated 
above -3.2 eV. These As derived bands are hybridized 
with the Fe d states. The Fe d states account for the 
bands between -2.2 eV and 2 eV, with La derived states 
occurring at higher energy. This is the ordering expected 
from the Pauling electronegativities: 3.44 for O, 2.18 for 
As, 1.83 for Fe and 1.10 for La. While there is some 
hybridization between Fe and As, it is not strong and is 
comparable to oxides. Thus a separation of the structure 
into independent LaO and FeAs subunits is not justified 
from the point of view of electronic structure or bonding 
and so we write the chemical formulae in order of elec- 
tronegativity. This ionic view is supported the similarity 
of the Fe d derived bands of LaFePO and LaFeAsO in 
spite of the chemical differences between P and As. This 
suggests that doping mechanisms other than replacement 
of O by F may be effective, i.e. on the La site (e.g. by 
Th) or the As site (e.g. by Te or Se). This may allow 
stabilization of the compound with a wider range of dop- 
ing levels than can be achieved using the O site, perhaps 
leading to higher T c . 

Fe d states would normally split into a lower lying e g 
manifold and higher lying t<i g states in an As tetrahe- 
dral crystal field. The gap between these would be at 
a d electron count of 4 per Fe. However, this crystal 
field competes with the direct Fe-Fe interaction to yield 
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FIG. 1: (color online) LDA density of states and projections 
onto the LAPW spheres on a per formula unit both spins 
basis. Note that much of the As p character will be outside 
the As sphere, reducing their apparent weight. 




FIG. 3: (color online) LDA Fermi surface of 
LaFeAsO shaded by velocity (blue is low velocity). 
The symmetry points are r=(0,0,0), Z=(0,0,l/2), 
X=(l/2,0,0),i2=(l/2,0,l/2),M=(l/2,l/2,0),A=(l/2,l/2,l/2). 
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FIG. 2: (color online) Band structure of LaFeAsO around 
Ef showing the effect As breathing along z by (*)2:as=0.04 
(0.035A). The unshifted band structure is indicated by the 
solid black line, while the shift away (towards) the Fe is indi- 
cated by the blue dotted (green dashed) lines. 



a more complicated band structure, with no clear gap at 
4 electrons per Fe. Rather, the main feature is a pseudo- 
gap at an electron count of 6. The Fermi energy for d 6 
Fe 2+ is at the pseudogap. While we do find sensitivity of 
the bands near Ep to the As height as shown in Fig. O 
the PBE band structure with the PBE relaxed structure 
is very similar to the LDA band structure, including the 
details of the Fermi surface. 

The Fermi surface (Fig. [3]). has five sheets: two high 
velocity electron cylinders around the zone edge M-A 
line, two lower velocity hole cylinders around the zone 
center, and an additional heavy 3-D hole pocket, which 
intersects and anti-crosses with the hole cylinders, and is 
centered at Z. The heavy 3-D pocket is derived from Fe 
d z states, which hybridize sufficiently with As p and La 



orbit als to yield a 3D pocket. The remaining sheets of 
Fermi surface are nearly 2D. The electron cylinders are 
associated with in-plane Fe d orbitals and have higher 
velocity and will make the larger contribution to the in- 
plane electrical conductivity, The Seebeck coefficients are 
positive because of the proximity to band edges in the 
hole bands: S xx =S yy =6.8 /iV/K and S zz =8.5 /iV/K, at 
300K. However, for N(Ep) and other quantities that de- 
pend on the density of states, such as spin- fluctuations, 
and electron-phonon coupling, the heavier hole pockets 
may be more important. Specifically, the three hole 
sheets together contribute 80% of N(Ep) but only 31% 
of N(Ef)v x ] . The average Fermi velocities are 0.81xl0 7 
cm/s (in-plane) and 0.34xl0 7 cm/s (c-axis) for the hole 
sections and 2.39xl0 7 cm/s (in-plane) and 0.35xl0 7 cm/s 
(c-axis) for the electron sections. Including all sheets, 
v xx =Vyy=1.30xl0 7 cm/s and u^=0.34xl0 7 cm/s. This 
yields a resistivity anisotropy of ~15 for isotropic scat- 
tering. We note that the DOS is rapidly changing near 
Ef and therefore these quantities will be quite sensitive 
to the electron filling, structure and other details. 

The volume enclosed by the two electron cylinders 
(equal to that enclosed by the hole sections) corresponds 
to 0.26 electrons per cell (0.13 per formula unit). Ef lies 
just above a peak in the DOS, which leads to a rapidly 
decreasing DOS with energy. This peak is associated 
with a van Hove singularity from the 3D hole pocket, 
which becomes cylindrical as Ef is lowered. The calcu- 
lated value at Ef is N(Ef)=2.62 eV _1 per formula unit 
both spins. The corresponding bare susceptibility and 
specific heat coefficient are xo=8-5xlO -5 emu/mol and 
7o=6.5 mJ/mol K 2 Thus LaFeAsO is a low carrier con- 
centration, high density of states superconductor. This is 
in contrast to the cuprates, which have high carrier con- 
centration (near half filling with large Fermi surfaces) 12 
and lower density of states. Recent experimental data 
also indicates low carrier concentration ] 13 i 14 i 15 Electron 
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(meV ^ FIG. 5: Phonon dispersions and density of states of LaFeAsO. 



FIG. 4: (color online) Induced moments as a function of an- 
tiferromagnetic (AF) and ferromagnetic (F) field. The field 
is applied by an extra spin dependent potential +/- B inside 
the Fe/Ni LAPW spheres. 



doping mainly shrinks the hole pockets, especially the 
3D sheet, since these have higher mass than the electron 
surfaces. This leads to an increasingly 2D Fermi surface. 
Within Stoner theory a ferromagnetic instability occurs 
when NI > 1, with TV now on a per spin basis. Since 
the DOS near Ep is close to pure Fe d in character it is 
appropriate to choose I ~ 0.7 - 0.8 eV, which would put 
LaFeAsO very close to a magnetic instability. 

Fixed spin moment calculations were done for the 
LaFeAsO, LaFePO and LaNiPO, and for LaFeAsO with 
O replaced by virtual atoms, Z=7.9 and Z=S.l (to simu- 
late doping) il& In the LDA LaFeAsO is indeed on the bor- 
derline of a ferromagnetic instability and electron doping 
moves away from this instability. The ground state is 
an itinerant ferromagnet with a moment of 0.08 /ig/Fe. 
However, the small energies involved are below the pre- 
cision of the calculation, so it may only be concluded 
that the material is on the borderline of ferromagnetism. 
This borderline behavior, where the energy is indepen- 
dent of magnetization, persists up to ~ 0.2 /i#/Fe. This 
is close to the value where the holes become fully polar- 
ized. This places stoichiometric LaFeAsO near a ferro- 
magnetic quantum critical point. 

The experimental susceptibility, 4 x(^) °f undoped 
LaFeAsO is weakly temperature dependent between 20K 
and 300K, with value ~ 50xl0 -5 emu/mol, and increases 
strongly at lower T. While the low T %(0) is not known, 
even taking the higher temperature value, the implied 
Stoner renormalization (1 — NI) -1 from Xexp/xo is 6 
and is probably significantly higher depending on the ob- 
served upturn below 20K. By comparison Xexp/x is ^5 
for MgCNi 3 ^ -7 for Sr 2 Ru0 4 A 8 and -9 for Pd metal. 
Ferromagnetic spin fluctuations are strongly pair break- 
ing for singlet superconductivity both in the s channel 
favored by electron-phonon interactions, and in d chan- 
nels, such as in the high-T c cuprates. Thus such fluctua- 
tions would suppress T c . However, it is conceivable that 
a sufficiently strong electron phonon interaction, which 



is needed to explain the T c ~ 26K in that scenario, could 
overcome the pair breaking effects of ferromagnetic fluc- 
tuations as was discussed for MgCNi3 J 17 i 19 Turning to the 
F doped material, which is the actual superconducting 
phase, Xexp is higher than for the stoichiometric material 
and increases as T is reduced over the whole tempera- 
ture range, rising to ~ 200xl0 -5 emu/mol at T c . This 
is opposite to the trends with band filling in N(Ep) and 
the fixed spin moment curves. The most simple expla- 
nation would be in terms of secondary phases. Another 
explanation is that the undoped compound is magneti- 
cally ordered. There is in fact a resistivity peak at ~ 
150K and a minimum at ~ 100K, 4 but x(T) does not 
show strong changes at these temperatures. 

We did calculations with an antiferromagnetic arrange- 
ment of the Fe spins in the square lattice. As for fer- 
romagnetic alignment, we find a borderline instability, 
with stable moments in the virtual crystal calculations 
at lower electron count while the proximity to antiferro- 
magnetism is reduced with increasing electron count. We 
also did calculations with small fields applied to the Fe 
LAPW spheres in both ferromagnetic and antiferromag- 
netic configurations. The resulting moments, defined by 
the magnitude of the spin polarization in each Fe LAPW 
sphere, are shown in Fig. [H along with results of similar 
calculations for LaNiPO and LaFePO. 

While ferromagnetism is initially favored, the antifer- 
romagnetic response becomes larger for higher fields with 
a cross-over at an energy scale of 30K. Thus thermal fluc- 
tuations at higher temperature and perhaps quantum 
fluctuations would both favor antiferromagnetic fluctu- 
ations over ferromagnetism. This nearest neighbor an- 
tiferromagnetism is related to super exchanged Doping 
with electrons rapidly reduces the proximity to ferromag- 
netism. Furthermore, it should be emphasized that La- 
NiPO, which is also a superconductor has much lower 
susceptibility for both ferromagnetic and antiferromag- 
netic fluctuations. 

The calculated phonon dispersions and density of 
states of LaFeAsO are shown in Fig. [5l The acous- 
tic modes are not strongly anisotropic. In particular, 
there are no soft elastic constants associated with shear- 



ing of the planes along the c-axis. The Debye tempera- 
ture from the acoustic mode velocities is 0£>=34OK with 
an uncertainty of 10% due to sampling. Turning to the 
optic modes, the dispersions may be divided into two re- 
gions. Below ~ 300 cm -1 the dispersions are dominated 
by metal and As modes of mixed character, while the 
modes above ~ 300 cm -1 are strongly O derived. The 
phonon density of states shows three main peaks below 
300 cm -1 . These are all of mixed metal and As charac- 
ter, and are a large peak centered at ~ 100 cm -1 and 
smaller peaks at ~ 170 cm -1 and ~ 280 cm -1 . 

The phosphides LaFePO and LaNiPO also supercon- 
duct, though with lower T c . The band structures near 
Ep of LaFeAsO and LaFePO (Ref. 11) are very similar. 
We find that LaFePO, like LaFeAsO is on the border- 
line of a ferromagnetic state in the LSDA. As expected 
from the different electron count, the electronic struc- 
ture of LaNiPO is very different. Calculations yield large 
two dimensional Fermi surfaces and lower N(Ep). With 
the assumption that the origin of superconductivity is 
the same in the Fe and Ni compounds, spin fluctuation 
mechanisms seem unlikely since these depend on a match 
between the q dependent susceptibility and the Fermi 
surface, requiring different coincidences for the Fe and 
Ni compounds. One scenario would be that there are 
strong electron phonon interactions related to the layered 
crystal structures and the presence of high charge ionic 
species in proximity to metallic layers. Then ferromag- 
netic spin fluctuations are pair breaking and the main dif- 
ference between the P and As compounds is the heavier 



mass of As, yielding a trend opposite to the usual isotope 
effect. Conversly, if the superconductivity of the Ni com- 
pound has a different origin (note the low T c of LaNiPO) 
a common spin-fluctuation based mechanism for the two 
Fe compounds is likely. Based on the suppression in the 
ferromagnetic susceptibility on doping the relevant fluc- 
tuations would be antiferromagnetic in nature. In both 
scenarios, doping plays a key role in the Fe based com- 
pounds. Electron doping reduces N(Ep), which would 
lower pairing strength by reducing the phase space, and 
at the same time strongly suppresses ferromagnetic fluc- 
tuations, which are strongly pair breaking for singlet 
superconductivity. Furthermore, antiferromagnetic spin 
fluctuations are also weakened by electron doping, which 
is of relevance if these are involved in the pairing. One 
way of distinguishing these pictures would be via dop- 
ing on the Fe site, e.g. with Co or Ni. Since the elec- 
tronic structure near Ep is from Fe d bands, such doping 
would be much more strongly scattering than doping in 
the O layer and would be very detrimental to non-5 wave 
spin fluctuation mediated superconductivity but not to 
s-wave superconductivity. It will be interesting to map 
out the similarities and differences between LaFeAsO and 
other superconductors with ionic metal type electronic 
structures, and systematically explore the relationship 
between spin-fluctuations and superconductivity. 
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